Hidden symmetry and Collective behavior 
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We study the relationship between the partially synchronous state and the coupling structure 
in general dynamical systems. Our results show that, on the contrary to the widely accepted 
concept, topological symmetry in a coupling structure is the sufficient condition but not the necessary 
condition. Furthermore, we find the necessary and sufficient condition for the existence of the partial 
synchronization and develop a method to obtain all of the existing partially synchronous solutions 
for all nonspecific dynamics from a very large number of possible candidates. 
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Synchronization has attracted extensive attention in 
the physical, biological, ecological, and other systems Q. 
The theory of synchronization focuses on the dynamical 
behavior of many-body coupled systems. Of the types of 
synchronization, global synchronization (GS) is a well re- 
searched topic Q. Partial synchronization (PaS), which 
is the synchronization only emerges in part of a system, 
is a more general synchronous phenomenon. And PaS 
has been observed in the systems whose parameters are 
outside the GS regime. A "simple" case, the PaS in 
a globally coupled system, is studied in detail in Refs. 
[J, |6[ . And random and other complex coupling struc- 
tures are studied in Refs. p-0]- However, the underlying 
mechanism of PaS is yet far from clear. A fundamen- 
tal difficulty lies in finding the PaS solutions for a given 
structure. And the problem remains open. 

Considering the dynamics [1(| 
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FIG. 1: (a) A topological structure of the coupled system 
without any partial synchronous solution, (b) The average 
distance di,2 and (£2,3 versus the coupling strength e with the 
structure shown in the inset graph. 



matrix Aa = 



X n+1 =F(X n ) + e(C<g>r)F(X r , 



(1) 




Here, nodes 1 and 2 are 



where X = (x 1 , x 2 , • • • , x w ) ' represents the state of a sys- 
tem consisted of N subsystems {x 1 }^ and "()'" means 
the transpose of a matrix. The independent state, x 1 7^ 
x 2 ^ • • • ^ x N , and the GS state, x 1 



N 



universally exist for common couplings |ll| . The PaS 
solution, however, may not exist in a given coupled sys- 
tem. As an example, for the structure shown in Fig. 
0Ja), all 203 candidate PaS solutions (e.g., x 1 = x 2 and 
x 3 ^ x 4 7^ x 5 7^ x 6 ; this is a kind of exhaust algorithm) 
do not satisfy Eq. (JTJ) . 

In recent years, it has been developed to a common be- 
lief that there exist close relation between the topological 
symmetry of a coupling structure and the PaS state. For 
example, an asymmetric PaS pattern that does not fol- 
low a symmetrical structure has never been observed [5j , 
the theory of symmetric group can be used to describe 
the periodic PaS state in several regular structures with 
the same symmetry Q , and all PaS states corresponding 
to each topological symmetry in a ring have been ob- 
served [3]. Based on these examples, one could suppose 
that symmetry is the necessary and sufficient condition 
for the existence of the PaS state. 

When the number of subsystems is small, the above 
statement may seem to hold true. The coupling struc- 
ture in the inset of Fig. [IJb) given by the adjacency 



/° 

1 

1 

Vi 

symmetric, and A4 is invariant under the permutations 
1 — > 2 and 2 — > 1. The curves d\p{e) and G?2,3(e) Q3 are 
shown in Fig. [T] (b), where e is the coupling strength. 
The synchronous solution x 1 = x 2 7^ x 3 7^ x 4 is observed 
in the region e e [0.3, 0.45] U [0.9, 1]. Thus, the PaS state 
will be achieved with the corresponding symmetry in A4 
among the symmetrical nodes. 

A more complex case is shown in Fig. [2] (a) with the 
same dynamics as above. There are two clusters and their 
nodes are denoted 1, 2, • • • , m and rii+1, ni+2, • • • , m + 
n-i- Node i (1 < i < n{) is coupled to nodes i — k,i — 
k + 1, • • • ,i + k, ni + i — l,ni + i — I + 1, • • • , m + i + /; 
and node rii + i (1 < i < ni) is coupled to the nodes 
ni + i — k, ni + i — k + 1, • • • , m + i + k, i — l,i — I + 
1, • • • , i + I. Obviously, there is "rotational" symmetry 
in every cluster; in other words, the adjacency matrix is 
invariant under a "rotation" permutation in each cluster; 
that is, 1 — > 2, 2 — > 3, • • • , n\ — 1 —> ni, n\ — > 1, and 

m + 1 -¥ ni + 2, m + 2 -* m + 3, • ■ ■ , 2m — l -> 2m, 

2m — s- n\ + 1. We define a M x M matrix Rm, where 
(Rm)m.i = 1, (#m)m+i = 1 (t = 1, 2, • • • , M - 1), and 
the other elements are 0. Thus, the permutation matrix 
[l3| of this transformation will be Td — R ni 0i?„ 2 , where 
© represents the direct sum of two matrices. Fig. [2] (b) 
shows the time-averaged variation in all subsystems <r(e), 
and that in the two clusters, u\{e) and 02(e) [12| for n\ = 
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FIG. 2: (a) The scheme of a topological structure with the 
"rotational" symmetry, (b) The variation of the two clusters 
o\, (72 and the whole system a as functions of e with m — 
n 2 = 100, k = 40, and I = 10. 



FIG. 3: (a) The scheme of an equal degree random struc- 
ture, (b) The variance of the two clusters ai, 02 and the 
whole network a as functions of the coupling strength e in 
the parameter n\ = 112 — 100, p — 1, and p r = 0.5. 



n 2 = 100, k = 40, and Z = 10. For e G [0.45, 1], <n = 
(T2 = 0; that is, Xi = • • • = x ni and x ni+ i = • • • = xat. 
The PaS solution of the "rotational" symmetrical nodes is 
observed. The above examples show that the presence of 
symmetry in a structure may be necessary and sufficient 
for the existence of PaS solutions. 

In this Letter, we investigate in detail the relation- 
ship between the PaS solution and the coupling structure. 
The PaS solution is denned as follows: For a dynamical 
system with phase space R 1 " 1 , a i^-cluster synchronous 
solution is a .RTm-dimensional subspace V of R Nm . It can 
be represented by 



X'l = X'l 

X J 2 = X* 2 



X'K 



X 'A- 



= X'l , 

i" 2 
= X l 2 , 

i n K 



(2) 



where x*« denotes the 6th subsystem in the ath clus- 
ter; and {rii\f =l are the sizes of each cluster, where 
Y,f= x ni = N (N > K > 1). The GS and indepen- 
dent solutions are the particular cases where K = 1 and 
K = N, respectively. The relationship between the PaS 
solution and the coupling structure could be described 
by the following two questions: 

Question A: If one finds symmetry in a coupling struc- 
ture, can a corresponding PaS solution be obtained? 

Eq. (|2|) can also be described by its matrix form: 



(T®I m )X = X,VXeV, 



(3) 



where T is a permutation matrix and I m is a m-dimension 
identity matrix. That is, X 6 V is invariant under the 
permutation transformation T ® I m , so V is the invariant 
subspace of T (g) I m and the eigen-subspace of T (g> I m 
corresponding to eigenvalue 1. That V is the invariant 
subspace of the dynamical system in Eq. ([1]) requires 



(c ® r)x g v, vx e v. 



(4) 



Next, if there is symmetry T in structure C, then C will 
be invariant under the permutation transformation T. 
The mathematical representation is 



or the matrices C and T are commutative in multiplica- 
tion. 

Question A can be represented by the following math- 
ematical statement: 

If Eq. © and Eq. © hold, then Eq. Q holds. 



© by X e V and 



T' X CT = C, 



(5) 



Multiplying the two sides of Eq. 
arbitrary T, we have 

[(TC) ® r]x = [(ct) ® r]x, vx e v. (6) 

Combining Eq. Q with Eq. (JSJ gives 

(T®/ m )(c®r)x=(C0r)x,vxev. (7) 

Thus, (C <S> r)X is also the eigenvector of T (?) I m with 
eigenvalue 1. Then, Eq. (g]) will be satisfied for all X G 
V. We conclude that for a symmetrical structure, the 
dynamical system has a corresponding PaS solution. 

Question B: If one finds a PaS solution, can the cor- 
responding symmetry in the coupling structure be ob- 
served? 

Here, an interesting example is shown in Fig. [3] (a). Let 
us consider two clusters; each has n subsystems and every 
subsystem is randomly connected to [pn] + 1 subsystems 
in the same cluster (p is a probability and [pn] means the 
integer part of pn) and [p r n] + 1 subsystems in the other 
cluster (p r is also a probability). All the subsystems are 
of equal degree and the connections between two clusters 
are also of equal degree. Fig. [3] (b) shows the variance 
of the two clusters (<7i and 02) and of the whole system 
(a) as functions of the coupling strength e for n = 100, 
p = 1, and p r = 0.5. As e G [0.34,0.76], a x =o 2 = 0, 
and a > 0, the PaS solution is observed. Due to the 
random connections between the subsystems, there is no 
symmetry in the structure. 

Question B can be represented by the following math- 
ematical statement: 

If Eq. © and Eq. g]) hold, then Eq. © holds. 

The following relations could be derived: 

(c <g) r)x e v, (r® i TO )x = x, vx g v, 

^ (r®/ m )(c®r)x = (c®r)x, 
<=* (T<s>i m )(C(s>r)x = (c®r)(T<s>i m )x, 
^ [(tc) ® (J m r)]x = [(ct) ® (r/ m )]x, 
^=^ [(TC - CT) ® r]x = 0. (8) 



Obviously, Eq. @ is not equivalent to Eq. ([5]). We can 
conclude that it is possible for a PaS solution to exist in a 
dynamical system without any symmetry; Fig. |3]presents 
an example. Another important point is that, in fact, the 
necessary and sufficient condition can be drawn from Eq. 
@ itself. 

For a simple representation but not loss its generality, 
we only discuss the case m = 1. The component form of 
Eq. US) is 

N 

Y, ( C kj - Qj )x» = 0, VX = (x 1 , x 2 , • • • , x N ) e V. (9) 
i=i 

for Tik = 1. By relabeling the subsystems to group sub- 
systems of the same cluster together, we rewrite Eq. @ 
as below: For cluster s, 



r N s +n„ — s 



y s , (a = l,2,— ,K), 

(10) 
where N s — when s = 1, and jV, = y^._ 1 rij when 
s = 2, • • ■ , /^ Q. Then Eq. © can be grouped into K 
terms as 

K N B +n s 

]T[ Y, (Cij-C kj )]y s = 0, (i = l,2,---,N). (11) 

s=l j=N s + l 

This is the ith row of Eq. ©. For all nonspecific dy- 
namics, we always have 

n s 

J2( C *.n s+1 - C KNs+] ) = 0, (s = 1, 2, • • • , *Q (12) 

forT it = l(i = l,2,.- ,JV). 

For a subsystem i which belongs to cluster s, we define 
the degree of subsystem i, which is contributed by the 
s'th cluster (s' ^ s), as the external degree of i from 
s'; And the degree of subsystem i, which is contributed 
by the sth cluster as the internal degree of i. Since 
Cij may be a fraction or zero, the external and internal 
degrees of a subsystem can also be fraction or zero. 

For subsystem i in cluster s, Y2j=i Ci.N a ,+j is the ex- 
ternal degrees of i from s' . And because of T^ = 1, both 
subsystems i and k are in cluster s. Then from Eq. (1121) . 
the external degrees of i and k from s' should be the 
same. And, Y^jLi Cj,JV s +j — Cm is the internal degrees of 
i. Because of Ckk = Cu = — 1> the rest part of Eq. (TT2"]) 
shows that the internal degrees of i and k should be the 
same. 

Considering the above two different situations in Eq. 
(1121) . we can have the following necessary and sufficient 
conditions of the existence of PaS solutions: SI: the 
external degrees of every subsystems in a cluster 
from another cluster should be the same. S2: the 
internal degrees of the subsystems in a cluster 
should also be the same. 

The two conditions form a complete representation of 
Eq. ([8]). Now we have a clear physical picture of the exis- 
tence of a PaS solution. And from the two conditions, we 



can derive all the PaS solutions in a given structure using 
the following strategy: The key is to divide the structure 
into several substructures, each comprising subsystems 
that can synchronize with each other but not with sub- 
systems in other substructures. Then the PaS solutions 
of the whole system will be combinations of all the pos- 
sible solutions of every substructure. 

To satisfy SI and S2, the PaS cannot be achieved be- 
tween directly connected subsystems whose degrees are 
coprime [15[ . So an N x N matrix S can be defined such 
that Si m — 1 if the degrees of the subsystems i and m 
are not coprime or d m — 0, and Si m = for other cases. 
Thus, the system can be divided into several groups of 
subsystems, and if the subsystems i and m are in different 
groups, then 5j m = 0. We can then suppose that x' — 
x m when Si m — 1- Thus, an equation like Eq.® can be 
represented as E J+e j + c '^ x ' 7+ + Ej_ e j_ C JJ^~ = °j 
where j + e J + = {j\Cf = C m j — Cij > 0} and 



3- G J- 



{J\C- 



(- a 



Cij < 0}. Positive (neg- 



ative) coefficients C^ (C*_) indicate that the internal 
or external degrees of subsystem m from subsystem j + 
(j'_) is greater (less) than that of the subsystem i from 
j + (]-)■ Therefore, according to the requirement of SI 
and S2, subsystem j + should belong to a cluster that 
includes one or more subsystems in J_ so that the posi- 
tive degree difference C + can be counteracted. But sub- 
system j + may not be synchronous with some members 
of J_ (some elements Sj + j_, Vj_ G J_ may be 0); if 

c t+ +T,j_eJ- S M- C J- > °: 3 J+ e J +> thcn Ec l-® c^- 
not be satisfied. This means x* = x m is impossible and 
Si m should be reset to 0. For the reasons specified above, 

if Cj_ +J2 3+eJ+ S J-J+ C t < °. 3 J- e J -> S "n should be 
reset to 0. To let i, m = 1, • ■ ■ , N and perform the op- 
erations discussed above, we obtain a new S. Thus, the 
system can be divided into some subgroups by consider- 
ing the new S. This procedure can be repeated until the 
new S equals the old S, and the size of the subgroups of 
the subsystems will decrease with every repetition. After 
repeated subdivision, we obtain a number of final groups, 
which can be grouped into clusters using the exhaust al- 
gorithm. If SI can be satisfied for every pair of clusters 
and S2 can be satisfied in each individual cluster, com- 
binations of these clusters are PaS solutions of the whole 
system. 

The number of candidate PaS solutions for any struc- 
ture rapidly increases with N |17l |. Thus, for a random 
structure with N = 1000 and p = 0.995 [H], it would 
be impossible to test all the candidates. But we can use 
the above procedure to find all PaS solutions. Three 
equal-degree groups that include at least two subsys- 
tems are found: group 1, {172, 261}; group 2, {532, 910}; 
and group 3, {231, 277, 503, 555, 756}. Thus, there are 
2 x 2 x 52 = 208 different PaS solutions. Fig. gja) shows 
the variances of the whole system, a, and the three clus- 
ters, <j\, CT2, and <T3, as functions of the coupling strength 
e. It is easy to find the PaS solutions xm — X2&\, 

^532 = X910, and X231 = £277 = ^503 = X555 = £ 756 . 
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FIG. 4: The variance of the PaS clusters and the whole sys- 
tem a as functions of the coupling strength e (a) A random 
structure: N = 1000, p = 0.995. (b) A Barabasi- Albert scale 
free structure: N = 1600. 



For the case of a Barabasi-Albert scale free structure 
(N = 1600) 0, Ell, Fig. EJb) show the existence 
of PaS clusters {102, 407}, {978, 1557}, {356, 815, 1035}, 
{266, 1142}, and {639, 1099} []]|. 

In conclusion, we studied the relationship between the 



coupling structure and the PaS state in general dynami- 
cal systems. The necessary and sufficient condition of the 
existence of a PaS state was found from an exact proof. 
And the results are counterintuitive; the existence of a 
PaS state does not require symmetry, as assumed pre- 
viously. And all of the candidate PaS solutions exist in 
a globally coupled system. Furthermore, as the exhaust 
algorithm cannot be used to obtain all of the existent 
PaS solutions, we developed a method to find all these 
solutions for a given structure. We focused mainly on 
the existence of PaS solutions in this letter. But to de- 
termine the stability of these solutions, the conditional 
Lyapunov exponent should be calculated by regarding 
the PaS manifold as the condition. Note that the proof 
also can also be applied to differential dynamical systems 
and the conclusions are the same. 
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The general form of T will be T = ®* 1 R ni (i = 

1, 2, • ■ ■ , K). Then the general form of X will be X = 

(y 1 ® li,m , y 2 <8> li.na , ■ ■ ■ , y K ® li,n K )' where 1 M ,N is 
an M x N matrix in which every element is 1. Then the 
Td in the second example is only a special case of this 
form, and the permutation transformation done to A4 
can also be represented by T4 = Rq © 1 © 1. 
Since dj in Eq.(|12[l can be a fraction, subsystems whose 
degrees are not coprime with each other may satisfy the 
requirement of this eqution (e.g., 1/6 + 1/6 + 1/6=1/4 + 
1/4). On the other hand, according to Eq. (|12p , a cluster 
with zero internal degree is also allowed. To satisfy the 
requirement of S2, subsystems in such a cluster cannot 
be connected with each other directly. Thus, two directly 
connected subsystems cannot be included in the same 
such cluster. 

The adjacency matrices of the two structures and a com- 
lete set of the PaS solutions of the second structres can 
be found in http://cluster.impcas.ac.cn/data.html 
For example, there are only 52 possible solutions (includ- 
ing the GS and the independent solutions) in a structure 
with five subsystems, but the number will be more than 
1.38 x 10 9 when the system size is just 15. 
A. -L. Barabasi, R. Albert, Science 286, 509 (1999). 



